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Abstract 

Eigenvalue correlations of random matrix ensembles as a function of an external 
(parametric) perturbation are investigated via the Dyson Brownian Motion model 
in the situation where the level density has a hard edge singularity. By solving a 
linearized hydrodynamical equation, a universal dependence of the density-density 
correlator on the external field is found. As an application we obtain a formula 
for the variance of linear statistics with the parametric dependence exhibited as a 
Laplace transform. 
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1 Introduction. 

Eigenvalues of N x N matrices can be viewed as energy levels {xa, a = 1, ■ • ■ , A^) of an 
effective Hamiltonian 

Ho = '^Xana, (1) 

a 

where Ha are the occupation at level a. A proposal by Wigner|]l| is that the levels are 
drawn from an ensemble of matrices, and when restricted to the eigenvalues has the joint 
probability distribution, 



P(xi, ■■■,XN)l[dXa = Cn e-'^^(^i'-'^-^") n dx„ (2) 



TV N 

n dxa = cn e-'^^(^i'-'^^'") n 

a=l a=l 

where 



W{xi, . . . , zat; m) = - ^ In |xa - + ^ u{xa), (3) 

a<h a 

Ctv is a constant and (3 = 2, 1,4 describes ensembles with unitary, orthogonal and sym- 
plectic symmetries. The density of levels is defined by 

— (E^(^-^«)) ' (4) 
\ a / eq 

where < ■ ■ ■ >eq is an average with weight P of Eq. (2). For u{x) = x"^ and Xa supported 
on the real line, a is given by the Wigner semi-circular law[|l| in the limit of large N , 

- x2, X G (-yiV,yiV). (5) 

Other ensembles arising from transport in disordered electronic systems 0, have Xa sup- 
ported on the right half line with different confining potentials, u{x). For u{x) = and 
a > 1/2 one finds the eigenvalue density 

/iV 

a{x) ~ W— , x^N (6) 



is universal near the "hard edge"0, in the sense that it is independent of a. It can be 
shown that for < a < 1/2, a (a;) ~ l/a;^^"|Q. In this paper we shall be interested in the 
response of the levels when is perturbed by an external potential. Of particular inter- 
est is the eigenvalue correlator as a function of the external potential. This problem was 
first studied in the context of the energy eigenvalues distribution of a disordered metallic 
ring subjected to an external magnetic field using diagrammatic techniques 0,0 and it was 
found that the eigenvalue correlations are universal after an appropriate rescaling. These 
results were later reproduced in0 using the Brownian motion model of Dyson[^ in the 
hydrodynamical approximation. Exact correlations for all strengths of the perturbation 



were obtained in^] using the method of supersymmetry pioneered by Efetov|[TO|l. All of 
the above results are valid in the bulk of the spectrum where the density is uniform; 
o"buik = c'"w(0) ~ constant. 

The phenomenological theory proposed by Dyson ^ interprets the eigenvalues Xa as po- 
sitions of classical particles which are governed by an over-damped Langevin equation 
subjected to Gaussian random force fjyj)^ 

dxa dW 

7-1- = -7^ + /a(r), (7) 
dr oxa 



where 7 is the friction coefficient, 



27 



fair) = 0, fa{T)MT') = -jSa^r ' t), (8) 

and T is related to the strength of the perturbation, X. Since the x^'s undergo a Brownian 



motion it is to be expected that X oc r P, O. A Fokker-Planck equation that describes 



^More precisely, the density of state-density of state correlation function at different fluxes. 
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the "time" dependent joint probability distribution can be derived, and reads 



dW ^ d 



P(xi,- ■ ■ ,xn,t), r > 0, 



(9) 



subjected to the initial condition; 

P{xi, - ■ ■ ,xn,0) = Y[S{xa 

a=l 

where x^ is the initial position of particle a. The stationary solution of the Fokker-Planck 



N 



equation is 



P{xi, ■ ■ ■ ,xn,oo) = Cn e 



It was shown by Dyson [H that the time dependent density, 



(10) 



(j(x,r) = -x„(r)) 



[here < • • • >^ denotes average with "time-dependent" weight given in Eq.(9)] satisfies a 
non-linear conservation law in the "hydrodynamical" approximation. The non-equilibrium 



density a{x, r) evolves in r according to 



d d ( d ^ 

cr(x, r) = — I a{x, t)^:^^! 



dr 



dx 



dx 



(12) 



where 



\E'(x, t) = u(x) 



dy(y{y, r) In \^-y\+{j^^ ^ \n[(y{x, r)]. 



(13) 



Note that the "time" dependent density is normed to A^; Jj^ dxa{x,T) = N, where K is 
the interval on which the levels are supported. 

The solution to Eq.(|l^) with Eq. (|T^) , will enable us to determine the parametric (r) 
dependence of quantities related to the eigenvalues Xa- 
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2 Linearization. 

The stationary density, cr{x) := a{x,oo), of the non-hnear diffusion equation satisfies a 
self- consistent Hiickel type equation, 



u{x) — J dya{y) In \x — y\ -\- ^— — ln[cr(a;)] = A = constant. 



(14) 



This suggests that for sufficiently long time|, we split the non-equilibrium density into an 
equilibrium part plus a small perturbation, p(x,r); 

a{x,T) = a{x) + p{x,t), (15) 

where a{x) the equilibrium density. Substituting Eq.(|15D into Eq.(p!2D and discarding all 
terms of 0{p^), gives 

d d 

— p(x, r) = —g^Ji^, (16) 

where the "particle" flux is 

J{x,t):= — / p{y, t) In \x -y\dy. (17) 

7 ox Jk 

This is unlike the ordinary diffusion equation in that the particle flux requires the entire 
distribution, p, to specify its value at one point. As an example, we consider the Gaussian 
ensembles with u{x) = x^, and K = {—^/N, VN). In the N ^ oo limit, and scaling into 
the bulk X <^ N, where the density is uniform, o"(x) = D = constant, the diffusion 
equation becomes, 

d , ^ D & 



— p(x,r) = - — — dyp{y,T)\n\x-y\. 



^More precisely Dt/^L ^ 1, where L is the interval over which the density extends | 



This is converted into 

^Pk{r) = \k\pk{T), (19) 

OT 7 

via a Fourier transform^, Pfc(r) = dxe^''^ p{x , r) , with the solution; 

Pkir) = Pk{r')e-^^'\'^^-^'\ r > r' . (20) 

Prom the r dependence of p^, we can infer the r dependence of 

Corr(x,y,r) := (a(x, r)(T(y, 0))eq - (a(x, r))gq (a(y, 0))eq (21) 
^ p^p^,-...-w-^|.|.Corr(p,fc), (22) 

where 

< <!> >eq:= (ll ^^°) ^^^i' ' ' ' ' " " " ' O 

denotes an average over the initial condition. Here Corr(p, k) is the Fourier transform of 
the equilibrium density-density correlation function|l|; Corr(x — y) = — , x ^ y. 

Q 

dx / cii/e*'=^+*^'^Corr(x - y) = —5{k + p)\k\. (23) 

-oo J-oo 7lp 



Therefore, 



Corr(x, 2/,r) = "5^3^- 



^ + tix-y) 



is the universal "time" dependent density-density correlation function when r is measured 
in units of ^ §, 0, 



The exact density-density correlation function for the Gaussian ensemble (/3 = 2, and x y) is 

2 



sin ■jr{x — y) 
TT{x-y) 



The hydrodynamical approximation essentially replaces [sin7r(a; — y)Y by 1/2. 
''[H gives the exact result for all t > and l3 — 2, 1, 4. 
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3 Hard edge correlations. 

We shall now focus our attention on ensembles where the eigenvalues are supported on 
the right half-line and (j{x) has the universal square root singularity at the origin. 

Corr(x, y, r) and lnt(a;, y, r) the twice integrated version of Corr(a;, y, r) are derived 

in Section 3. These are applied to determine the parametric dependence of the variance 
of arbitrary linear statistics, given by Eqs. (51) and (54). As an application, we compute 
in Section 4 the variance of J2a[^ + Xair)]^'^ , which gives the conductance fluctuation 
of a quasi-one dimensional disordered system as a function of the external perturbation. 
According to the Landauer formula, the conductance, g, is J2a[^ + ^a]~^- 
We have therefore solved in the hydrodynamical approximation the problem posed in[^ for 
the case where the equihbrium eigenvalue density displays a hard edge singularity; a{x) ~ 
D I y/x. Since translational invariance is no longer valid, results obtained in the bulk scaling 
limit of the Gaussian ensembles (with u{x) = x^) are no longer applicable here |T^] . For 
example, the gap formation probability, -E^(0, s), of the Laguerre ensembles (with u{x) = 
X — alnx, X > 0, a > —1.) is distinct from that of the Gaussian ensembles[^, |14 . 



We expect Corr(a;, y, r) to have distinct "time" decay modes from that found in the bulk 
scaling case 0. 

^An application of a dimensional argument on Eq. (18), obtained by scaling into the bulk of the 
spectrum, shows that the typical distance covered by a diffusing particle over time t is \x\ ^ t, which 
is "faster" then Einstein diffusion; \x\ ^ \/t. Looking ahead to Eq. (24) the same analysis shows that 
|a;| ^ t^^^. This suggests that particle transport near the hard edge is intermediate between ballistic 
motion and classical diffusion. 
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The diffusion equation now reads 



dt dx \ ^Jx Jo X — y 

where Dt/'j = t. [We have written p(x, r) = p(x, t).] With the ansatz p{x, t) oc p(x, A)e~'*'*, 
Eq.(24) becomes, 

Ap(x.A) = #ardy^V A>0. (25) 



dx \ Jo X — y 

The boundary condition on p{x, A) is such that the particle flux vanishes at the boundaries; 
lim^_,o J{x, A) = 0, and limj-^oo J{x, A) = 0. The boundary condition at x = reads, 

lim^Pr^ci, = 0. (26) 
x^o y/x Jo X — y 



Since P y/y(x-y) ~ ^i 2; > 0, we find it convenient to write p(x. A) as 



p(a;,A)=p(x,A) + ^, (27) 

'X 



where p(x. A) fulfills a stronger condition, 

r,,*^ = 0. (28) 

Jo y 

and est (A) to be determined later is a function of A only. In appendix ^ it is shown 
that the final solution (Eq. (42)) fulfills the boundary condition (^). The term ^^^^^ 
is similar to the equilibrium solution and we conjecture that the solution for a general 
hard edge density always has this structure. We now go on to solve Eq. (EH). This can 



be accomplished by "un-folding" the half-line into the real-line |T5|. With the change of 
variables, v = ^/y, u = y/x, we have 

~ iyEi^ ^ r , pijtR = f A, (29) 

X — y Jo X — y Jo w' — v' 



dv "^^ ' ' + / dv ' \ (30) 

u — V Jo U + V 



Now introduce an odd function of v 



and find 



pi(^7,A) = P{v\X), v>Q (31) 
Pi(t;,A) = -p(^;^A), < (32) 



X — y Jo u — V Jo u + V 

- d,MEi^ + f <i.„^ = r d„^i^. (33) 



U — V J-oo U — V J-oo U — V 

With Eq. (|33|) and Eq. (p7|) the diffusion equation, Eq. (pSj), becomes, 

cst(A) 19/1 /•+°^ J Pi(t;,A)\ 

+ Api(M, A) = — — - / dv "^'' . (34) 



u 2udu \uJ~oo u — v^ 

Although Eq. (p4D is derived for m > it is also valid for u <0. This can be seen from the 
fact that pi is an odd function of u. To proceed further, we introduce the even function, 

P2(m,A):= , (35) 

u 

and with its aid we find p2(^, A) satisfies, 

2Acst(A) + 2Au2p2(M,A) = ^ (35) 

OU J-oo U — V 

We have made use of the flux condition at origin, Eq. (|28|) , 

30 POO 

dup2{u, X) = 2 / dup2{u, X) = 0, (37) 







to arrive at Eq. (|36|) . This accomplishes the un- folding of half-line diffusion equation onto 

the real line. Clearly, Eq. (|36|) can be be solved by a Fourier transformation. 

Eq. (|37D implies that P2{u, A) is an oscillatory function of u and its Fourier transform, 

P2(A;, A) := due^^'^p^iu, A), (38) 



vanishes at /c = 0. A further condition is that p2 is even in k. 

Once p2{k,X) is found, the original density can be recovered by standard inversion for- 
mulas; 

p{x, A) 1 



J dk cos (k\^^ p2{k, X). (39) 



X TC Jo 

A simple calculation shows that the transformed density, p2{k, \), satisfies an Airy equa- 
tion with a point source, 

^p,{k,X) + ^P2{k,X) = 27rcst{X)5{k). (40) 

The solution of this is a linear combination of Ai (7r/2A)^/^|/c|^ and Bi ^— (7r/2A)^/^|A;| 
and reads up to a constant. 



P2ik,X) = 3^/6r(2/3) (2A7r5^'^' 



Bi(0)Ai ( - ( ^V^' |A;| ) - Ai(0)Bi ( - f^V^' I'^l 



2A/ ' 7 V V2A 



^vi*i^./.(|(^)"Vn. (41) 

where we have made use of Ai(0) = ^2/3) ' -^KO) = ^"^^ relations between 

Ai(— a;), Bi(— s) and J±i/3 (^^x^^"^^ to arrive at Eq. (p]). Note that p2{k,X) vanishes 
at = and is even in k. We now use the jump discontinuity of ■^p2 across = to 
determine est (A): 



3-2/3 Ta:/n^ _ 3-1/6 



c=t(A) = - 



fc=-0 



TT \ 6 



18 A 



Transforming back to p{x, A) and with the help of Eqs. (|39| ) and (^Tj) we get 

p(x,A) = r dkcosiky/^)Vk3i(l\f^\k\l) + ^^^. (42) 
JO 3\3V2A / ^yx 

Therefore the "time" dependent density reads, 

PCX) 

p{x,t)= I dXC{X)e-^'p{x,X), (43) 
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where C(A) is determined by an initial condition. 

Eqs. (42) and (43) will now be used to find the dynamical density-density correlation 
function. It is clear that Corr(x, y, t) and 

ry 

lCoTi{x,y,t) := / dz Cott{x, z,t) (44) 
Jo 

satisfy the same diffusion equation as p{x,t). Hence, 

roc 

lCoTT(x,y,t) = dXC(X,y)e-^' p(x,X), (45) 
Jo 

and 

PX 

lnt(x,y,t) := / dz lCorT(z,y,t) (46) 
Jo 

= -kJ^ dXC{X,y)e~^' dk Vk sm{k,/^) JyJ- i^—j p/2j,(47) 
where we have used formula ( [7^ ) in appendix ^ to arrive at the last equation. C{X,y) 



is determined by the "initial condition", lnt{x,y,t = 0) = ^ In 



\Vx-Vy\ 



A derivation 



of this formula can be found in appendix 0. To find C{X,y) we first perform a Fourier 



sine transform in the variable ^/x^, /q°° d^/x sm{py^) ■ ■ ■ , on Eq. (^) at t = to get an 
integral equation satisfied by C{X,y)f\ 

-^-(pVi?) = yi imx,y)^v.h„{j^/''') . (48) 

This can be rewritten as 

2sin(p^) _^ (C{i^.,y) 2_ 



^1/3 ^fV^, If" , (49) 



/57rV/' I u"" '3 



^We use the transform Jq°° dxlnjl^ll sin(px) = ^ sin(pa). 
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where T-Cu{f{u),^) = /^(O •= Sq^ duf{u)uJ,j{u^), denotes the Hankel transform. Using 
the Hankel inversion theorem [/(x) = Jq° dC^fuiO^J'^i^'^)^ see[jl6[], we find, 



^ dpsm{p^)y/pJi/s (^gP^^^w) • (50) 



3/37r4 Jo 

Therefore the "time" dependence of the twice integrated correlation function is displayed 
as a Laplace transform with the spectral parameter A := 

Int(x, y, t) = --— / -e-^* F{x, \)F{y, A), (51) 
Spn^ Jo A^ 

where F{x, A) := Jq dzp{z, A) and p{z, A) given by Eqs. (27) and (42). Thus the variance 
of an arbitrary linear statistics, 

Q{t) ■■= J2QMt))= dxJ26{x-Xa{t))Q{x), (52) 



IS 



Var[Q,t] := < Q(t)Q(0) >eq - < Q(t) >eq< Q(0) >eq (53) 



oo roo roo roo 

dx dyQ{x)Q{y)CoTT{x,y,t) = dx dyQ'{x)Q'{y)lnt{x,y,t). (54) 
Jo Jo Jo 



This is the main result of our paper. 



4 Variance of -ri—. 

l+x 



As an example we take Q{x) = which is the Landauer formula for the conductance. 
A simple calculation gives 



Var 



1 

t 



with 



Jo 1 + xr 
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1 

2 JO 



5/6 oo ^18A\"^^ 



n=0 



TT 



(56) 



where 



r(l + n/3) 



r(-(n- l)/3)' 



(57) 



Eqs. (56) and (57) can be found in|[Tj]. Now when n = 1 mod 3, An is taken to be zero. 
The remaining terms can be grouped into two sums with n = mod 3 and n = 2 mod 3 
respectively. An apphcation of the ratio test to these series shows that both converge 
absolutely and uniformly for all A. We deduce that (with j3 = 2), 



^ n=0 



(5^ 



converges uniformly in the A^^'^ plane. Using this representation and after performing the 



Laplace transform, we get, 
1 



Var 



where 



1 + x 



1 oo Y ( ndl\ 

^2/3 



n=0 



t3 



<-^2/3 _|_ <-^4/3 _|_ <-^5/3 _|_ *-"2 _|_ <-^7/3 ^ 



^2/3 ^4/3 ^5/3 ^2 



^7/3 



(59) 



(^)'^'r(2/3) 

C2/3 = ^ L.,. , J ~ 0.050346, 



2r2(i/3) 



3(^)'^'r(5/3)r(4/3) 

r(-i/3)r(i/3) 



a 



4/3 



C5/3 = 
C2 



3f^)^/'r(5/3) 



r(i/3)r(-2/3) 

27r2(5/3) 



-0.0632871, 



0.042735, 



7/3 



-r(7/3) 



47r2r2(-l/3) 
9(^)^/'r(5/3) 9^^ 



0.0337737, 

1/3 



47r7 



r(8/3) 



r(-i/3)r(-2/3) ior(-4/3)r(i/3) 



-0.0716263. 
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This completes the determination of the long "time" behaviour of Var[l/(1 + x),t]. Ac- 
cording to the Watson Lemma [p!^, Eq. (59) is an asymptotic expansion. The fact the 
Brownian motion ensemble tends to the stationary ensemble with exponential speed as 
soon as tD/jL ^ 1, enables us to conclude that limf^oo Var[Q, t] = 0. At t = 0, we simply 
go back to Eq.(55) and make use of the completeness theorem of the Hankel transform to 
show that Var[l/(1 + x), 0] = a result obtained previously 0, |20|, Other examples 



of Q can be found in! 19 
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A Check of solution 



In this appendix we are going to show that 



p(x,X) = v^cos(A;v^)v^Ji f-^/^|A;|tVA; (60) 

\/x Jo 3 V 3 V 2 A / 



is a solution to the integro-differential equation 

To obtain better convergence in the manipulations we will instead show p{x, A) fulfills 
the equation obtained by integrating Eq.(|6T|) from to z. From the boundary condition 
Eq. (^) we see that the lower boundary term on the right hand side to vanish. Using this, 
Eq. (|61|) becomes 

\ r dxp{x,\) = r ^^^dy. (62) 

Jo y/Z Jo z — y 



This equation is a stronger statement than Eq. (|6T|) since any constant will vanish when 
Eq. (|62D is differentiated to give Eq. (0). We start by evaluating the right hand side of 
Eq.(62), 

R.H.S = ±r P^dy (63) 
^z Jo z-y 

d.ii 

+ 



1 r (ifA)' dy 



, 1 r r ^^^(^ v^) Vkj^ f 2 g |,|f ) dk 

^JzJo z — y y/y ^/z Jo Jo z — y a\3V2A / 



Using the identity, 

p ^cos{ky) ^ , r ■ r\ 

/ dy = ZTT oik) + TT ^/zsm[k ^/Z), 

Jo z — y 
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in Eq. (64) gives 



1 f°^p{y,X) 



z Jo z — y 



dy 



sm{k^)Vk3. {^-^ dk. (65) 



We now turn to the left hand side of Eq. 



L.H.S 



A / dxp{x,X) 



(i8a) 



(66) 



2A 



'X 

1 

TT \ 6 



+ A / dx [ \fx cos{k y/x) \Pk Ji {-S—r \k\'^\dk 
Jo ^\3V2A 



TC 



\pz + 2 A / / dtt^ cos (A;t)v^Ji -Jtt- |A;|n rfA;.(67) 
18 A/ Jo Jo 3 V3 V 2 A 



The integral with respect to t = y/x is 



dtt^ cos{kt) 



d'^ j sin(/c yi) 



dk^ 



k 



The left hand side is then 



L.H.S 



2A 
+2 A 



v^-2A 



18 A 

°° d j sm{k y/z) 1 d ( 
dk\ k j dk 



■2 nr 

t V3 V 2A 



dk. 



(69) 



The boundary terms from the partial integration are zero, since 



VkJi f Q \/ 77T l^h ) ~ const /c 4 cos(const/c2 — phase), k oo 
3 V o V 2 A 



0, k = Q. 



(70) 



Now, 



L.H.S 



2A 



vr \ 6 



18 A 



v^+2A 



sin(fcv/i)|_r^ /2 , 
A; 9A; I 3 V3V2A ' ' 



-2 A 



sin (A; y^) 9' 




dk. 



(71) 
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The boundary term at infinity vanishes. The boundary term at the origin is 



k=0 



(72) 



By differentiation, 



2A HsV 2A ' ' 



(73) 



Gathering the pieces we have 



A 



/ dxp{x,X) = 71 [ sin{k ^/z) Vk Ji ('^\/'^ \k\2\ dk 
Jo Jo 3V3V2A/ 



A / dxp{x,X) 
Jo 



z Jo z — y 



-dy 



(74) 
(75) 



Here we make two remarks on the above derivations: 



1. Since the number of particles in the interval (0, z) goes to zero, when z tends to zero, 
the last equation shows that the zero flux condition at the origin is fulfilled. 

2. This is concern with the z —* 00 limit in Eq. ([7^) . The integrand on the right hand side 
is the product of two strongly oscillating functions, which oscillate out of phase. Since 
none of them are absolutely integrable, a Riemann-Lebesque lemma can not be used 
to conclude that the integral tends to zero. We may appeal to the theory of generalized 
functions that it is 0. The R. H. S. of Eq. (74) can be re-written up to irrelevant constants 



as 



J{a) = dy sm{ya)y/yJi/3{by^'^), 



where a := y/z and & := | (^) ^ • We wish to determine J'{a) in the limit a 
b{> 0) fixed. Consider instead. 



00, with 



y 



5 so 5 
4 2 
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The integral we need can be defined as the analytic continuation of J^{a) to ji — —1/2. 
Now, since lima_^.oo ^^^^^^^ = 5{x), we conclude by integrating over the 5 function that 
linia^oo = 0. Hence J{oo) — 0. Physically this means that the total number of 

particles in each lambda mode is zero. The zero mode, of course, contains particles. 
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B The equal time variance of a linear statistic. 

The equal time variance of a linear statistic Q is 

/•oo poo 

Var(Q,0) = dx dy Q{x) Q{y) CoTT{x,y,0), (76) 
Jo Jo 



a result first derived in by the method of functional derivative. The claim in is 



that the method is valid for all potentials, however, this is flawed by the use of the N = oo 
density, which may not exist for all potentials. For the sake of completeness, we present 



in this appendix a small modification of the arguments given in ||I9| and show that the 
result has a more general validity. Instead of focusing on Eq. ([76|), we prefer to work with 
the partially integrated expressions 

roo roc 

Var(Q,t) = - / dx dy Q{x) Q! {y)lCoii{x,y,t), (77) 
Jo Jo 

and 

roo roo 

Var(Q,t) = dx dy Q'ix) Q'iy)lntix,y,t), ' (78) 
Jo Jo 

where 

lCoTT{x,y,t) = / dzCoTiix, z,t), (79) 
Jo 

rx rx ry 

lnt{x,y,t) = / dz ICoii {z, y, t) = / dzi / dz2CoiT{zi, Z2,t). (80) 
Jo Jo Jo 

We will start by finding ICorr(x, y, 0). By definition crgas(a;) = J2'^=i — x°) and 

J,^dxl..J,^dx%ag^s{x)e~'^^-^ ^^^^ 



< crgas[x) >eq 



J,^dxl...J,°°dx%e-^l-] 



^To avoid the boundary terms it has been assumed that the hnear statistic fulfills -^/x Q{x)\^^q = 0. 
The linear statistic corresponding to the conductance satisfies this criterion. 
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The dependence of the energy, W, on the coordinates x?, ■ ■ ■,x% has been suppressed. 
The external potential u{x) is assumed to be bounded at the origin which is the physically 
interesting situation. Now write u{x) — — Jq f{z) dz + m(0), where f{x) is the force. The 
constant u{Q) can be set to zero as this can always be accomplished by a redefinition of 
the zero point energy. The energy is now 

W[u] = W - f{z)dz =-V/ f{z) eix'^^ - z)dz- y^ln\x^. -x% (82) 

where 9 is the Heaviside step function. The functional derivative of W[— /q f{z) dz] with 
respect to / is 

M^MiM . .gas(..0)d.. (83, 

The functional derivative of aeq(a;) =< (Jgas{x,0) >eq with respect to / is recognized 



as ICorr(x, y, 0), 



The density aeq{x) is approximated by 



1 b{NJ)-x rmJ) dz I z 
'^^"''^'^^ = -^^i X Jo V^i biNJ)-z ^^'^^ ^''^ 

where b{N,f) is the upper limit of support of a{x,N;f). ^'^g^/^^'^^ is computed as the 

linear term in e, when e5{x — y) is added to the force f{x). We have 



where r) — r]{€, y,N,f) — e ^ Jq ^ ^J^^dx — —e^ is the number of particles 
associated with the force e5[x — y) and h := h{N, f + eS{x — y)) = b{N — r], /) is the upper 
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limit of support of a{x; N, f + e6{x — y)) and a{x, N — rj; f) 0. The functional derivative 



IS 



6a{x,N) 



1 b — x^l y da{x,N) 



Ptt'^ \ X X — y b — y Ptt b — y 



dN 



In the limit oo, = V and ^'^jjjf^^ tend to zero. The extra density, Aa{x), 

when e particles is added to the system, experiences no force in the interval (0, b). There- 
fore Acr(x) is less than the corresponding density with e particles in the box (0,6): 
^— ^===^, X G (0,6). Outside (0,6), Acr{x) tends monotonically to zero. We conclude, in 
this hmit, the functional derivative is reduced to 



6a{x,N) 
ICorr(x, y, 0) 



1 fy 1 



(3n'^ \ X X — y 



' tin 



/Stt^ dx 



X 



(89) 
(90) 



From this it is seen that 



Int(a;, y, 0) 



In 



(91) 



and 



poo poo 

Var(Q,0) = dx dy Q' (x) Q' (y) \n 

pvr^ Jo Jo 



V^ + Vy 



(92) 



Thus at equal time, variances are independent of the potential u. 



They have the same upper hmit, because the part of (7{x,N;f + ed{x — y)) corresponding to /: 



ZJE dz 
c Jo z—x 



■ f{z) is a sohition, with N ~ rj particles, to the problem with the external force 



/, obeying the right boundary conditions. In other words — ^ \/ ~r' lo 'z^ \l /('^) ^ "^(^j ^ ~ ^5 /)• 
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